function eta = proj_tangent_space_stiefel(eta,U)
% This defines the projection
% of direction \eta in R^{n \times p} onto the tangent space of U
%
% Authors:
% Bamdev Mishra and Gilles Meyer
% {b.mishra,g.meyer}@ulg.ac.be

  eta = eta - U*symm(U'*eta);